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Abstract 

This  objective  of  this  work  was  to  predict  the  physicochemical  properties 
of  six  energetic  materials.  The  energetic  materials  of  interest  were  2,4- 
Dinitroanisole  (DNAN),  3,4-Dinitropyrazole  (DNP),  n-Methyl-p- 
nitroaniline  (MNA),  i-Methyl-2,4,5-trinitroimidazole  (MTNI),  3-Nitro- 
i,2,4-triazol-5-one  (NTO),  and  1, 3, 5-Triamino-2, 4, 6-trinitrobenzene 
(TATB).  The  physicochemical  properties  of  interest  were  the  aqueous 
solubility,  the  octanol-water  partition  coefficient  and  Henry's  constant. 
Three  reference  compounds  were  chosen  to  provide  validation  of  the 
proposed  approach:  cyclotrimethylene-trinitramine  (RDX),  octahydro- 
i,3,5,7-tetranitro-i,3,5,7-tetrazocine  (HMX),  and  8-CL20. 

Group  contribution  methods  available  in  the  literature  were  used  for 
predictions.  In  addition,  a  number  of  physical  properties  were  required, 
including  critical  temperature  and  pressure,  normal  boiling  point,  melting 
point,  heat  of  fusion,  and  vapor  pressure.  These  were  also  predicted  using 
the  chemical  structure  of  the  energetic  compounds.  Where  literature  data 
were  available,  comparisons  with  predictions  were  performed. 

Another  purely  predictive  method,  based  on  chemical  structure,  was  also 
employed  for  the  estimation  of  physicochemical  properties.  This  method, 
based  on  the  “Conductor-like  Screening  Model  for  Real  Solvents,”  required 
quantum  chemical  calculations  performed  using  TURBOMOLE  to  obtain 
the  optimized  electrostatic  potential  surface  of  each  energetic  material. 

The  resulting  surface  was  then  used  within  the  COSMOtherm  software  to 
predict  solution  behavior,  including  aqueous  solubility,  octanol-water 
partition  coefficient,  Henry's  constant,  and  vapor  pressure. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to 
be  construed  as  an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 
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1  Introduction 

Background 

The  use  of  chemical  structure  as  a  means  to  predict  physical  and  transport 
properties  and  physicochemical  properties  has  been  the  subject  of  intense 
research  in  recent  years.  Physical  and  transport  properties  examined  in¬ 
clude  melting  point  and  density  (Trohalaki  et  al.  2005),  and  viscosity,  sur¬ 
face  tension  and  thermal  conductivity  (Kauffman  and  Jurs  2001).  Physico¬ 
chemical  properties  examined  include  the  octanol-water  partition 
coefficient  (Marrero  and  Gani  2002;  Derawi  et  al.  2001),  aqueous  solubili¬ 
ty  (McElroy  and  Jurs  2001,  Ran  et  al.  2001;  Liu  and  So  2001),  soil- 
sorption  coefficient  (Delgado  et  al.  2003),  vapor  pressure  (Basak  and 
Mills,  2001),  and  Henry’s  law  constant  (English  and  Carroll  2001).  Two 
primary  approaches  have  been  taken  in  the  prediction  of  physicochemical 
parameters  based  on  the  structure  of  a  material:  (1)  group  contribution 
methods,  and  (2)  quantitative  structure-property/activity  relationships 
(QSPR/QSAR).  The  group  contribution  methods  are  primarily  based  on 
the  numbers  and  types  of  molecular  groups  that  form  the  compound  (Pol¬ 
ing  et  al.  2001);  thus,  the  first  order  group  contribution  methods  fail  to 
distinguish  between  isomers  that  contain  the  same  numbers  and  types  of 
groups.  Second  order  group  contribution  methods  incorporate  interac¬ 
tions  between  next  nearest  neighbors  in  an  attempt  to  address  the  exist¬ 
ence  of  isomers  and  their  influence  on  the  properties  of  interest. 

Quantitative  structure-property/activity  relationships  connect  physical  or 
chemical  properties  to  a  set  of  molecular  descriptors.  An  extensive  body  of 
research  in  this  area  has  focused  on  development  of  these  relationships  for 
use  in  the  pharmaceutical  industry,  where  aqueous  solubility  and 
lipophilicity  are  used  to  estimate  the  pharmacokinetics  of  potential  drugs 
(Horning  and  Klamt  2005).  QSPR/QSARs  also  have  potential  applications 
in  the  fuel  science  field  (Katritzky  and  Fara  2005). 

The  molecular  descriptors  can  be  constitutional  (depend  on  compound’s 
composition),  geometric  (depend  on  3D  orientation  of  atoms  in  compound), 
electrostatic  (depend  on  partial  charge),  or  quantum-chemical  (depend  on 
compound’s  electron  wave  functions)  (Katritzky  and  Fara  2005).  What  is 
critical  in  the  development  of  QSPR/QSARs  is  the  identification  of  the  ap¬ 
propriate  set  of  descriptors  that  allow  the  desired  attribute  of  the  compound 
to  be  adequately  predicted.  A  key  limitation  to  the  use  of  QSPR/QSARs  is 
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that  the  set  of  compounds  used  to  develop  the  relationship  should  be  similar 
to  those  compounds  for  which  predictions  are  desired. 

A  third  approach  (Klamt  and  Eckert  2000)  provides  for  a  priori  predictions 
of  thermophysical  data  for  chemical  systems.  Based  on  “Conductor-like 
Screening  Model  for  Real  Solvents,”  COSMO-RS  has  a  strong  quantum 
chemical  and  statistical  thermodynamics  basis,  which  allows  a  physically 
meaningful  description  of  the  interactions  between  molecules  in  a  solution. 
It  has  been  used  to  predict  a  number  of  physicochemical  properties  includ¬ 
ing  aqueous  solubility  (Klamt  et  al.  2002),  soil  sorption  coefficients  (Klamt 
et  al.  2004),  gas  solubilities,  and  partition  coefficients  (Maassen  et  al.  1995). 

These  methods  are  challenged  when  little  or  no  experimental  data  are 
available,  often  the  case  for  particularly  complex  compounds,  or  com¬ 
pounds  in  the  early  stages  of  development.  One  must  develop  a  means  to 
validate  the  predictions,  often  by  predicting  similar  compounds’  properties 
for  which  experimental  data  are  available.  Compounds  chosen  to  serve  this 
function  in  this  work  were  S-CL20  (2,4,6,8,io,i2-Hexanitro-2,4,5,8,io,i2- 
hexaazaisowurtzitane),  cyclotrimethylene-trinitramine  (RDX),  and 
octahydro-i,3,5,7-tetranitro-i,3,5,7-tetrazocine  (HMX).  All  three  struc¬ 
tures  contain  multiple  nitro  functional  groups  and  also  contain  cyclic  ni¬ 
trogen  groups.  The  structures  of  the  compounds  of  interest  in  this  work, 
which  contain  these  functional  groups  —either  cyclic  nitrogen  atoms  or 
nitro  groups  attached  to  the  ring,  or  both—  are: 

•  dinitroanisole  (DNAN) 

•  n-methyl-p-nitroaniline  (MNA) 

•  nitro-triazolene  (NTO) 

•  triaminotrinitrobenzene  (TATB) 

•  dinitro  pyrazole  (DNP) 

•  and  m-trinitro  imidazole  (MTNI). 

Objectives 

The  objective  of  this  work  was  to  develop  a  means  to  validate  the  predic¬ 
tions  of  physical  and  transport  properties,  and  of  physicochemical  proper¬ 
ties  through  prediction  of  similar  compounds’  properties  for  which  exper¬ 
imental  data  are  available,  specifically,  8-CL20  (2,4,6,8,10,12-Hexanitro- 
2,4,5,8,10,12-hexaazaisowurtzitane),  RDX,  and  HMX. 
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Approach 

This  work  was  accomplished  in  the  following  steps: 

1.  Six  compounds  of  interest  in  this  work  were  specified  and  three  reference 
compounds  along  with  their  Chemical  Abstracts  Service  (CAS)  registry 
numbers  were  identified. 

2.  The  molecular  structure  of  each  compound  was  used  to  generate  a  map¬ 
ping  of  atoms  in  each  molecule  in  a  Cartesian  coordinate  system. 

3.  Physical  properties  for  the  compounds  of  interest  were  estimated  and  the 
reference  compounds  was  undertaken.  In  addition  to  the  conventional 
physical  properties  estimated  from  molecular  structure,  the  Abraham  de¬ 
scriptors  were  computed. 

4.  Physicochemical  and  physical  properties  for  the  compounds  of  interest 
were  also  estimated. 

5.  Based  on  these  estimates,  a  variety  of  treatment  technologies  was  exam¬ 
ined  for  the  remediation  of  waste  streams  containing  energetic  com¬ 
pounds. 

Mode  of  technology  transfer 

This  report  will  be  made  accessible  through  the  World  Wide  Web  (WWW) 

at  URLs: 


h  ttp :  /  /www.  cecer.  a  rm  v.  m  i  I 
http://libweb.erdc.usace.armv.mil 
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2  Computational  Approach  and  Results 

Molecular  structures 

Table  l  lists  the  six  compounds  of  interest  in  this  work  and  three  reference 
compounds  along  with  their  CAS  registry  numbers.  The  reference  com¬ 
pounds  were  chosen  based  on  the  availability  of  literature  data  that  could 
be  used  to  validate  the  estimation  approach  and  their  inclusion  of  relevant 
molecular  groups,  in  particular,  the  nitro  (-N02)  molecular  group  (Figure 
l).  Figure  2  shows  the  molecular  structures  for  the  six  compounds  of  inter¬ 
est  and  the  three  reference  compounds,  obtained  from  the  Scifinder  Schol¬ 
ar  database  (ACS  2007). 


Table  1.  Compounds  of  interest  and  reference  compounds. 


Compound  Name 

Acronym 

Reference 

(Yes/No) 

Molecular 

Formula 

CAS  Registry 
Number 

2,4-Dinitroanisole 

DNAN 

C7H6N2O5 

119-27-7 

3,4-Dinitropyrazole 

DNP 

C3H2N4O4 

38858-92-3 

n-Methyl-p-nitroaniline 

MNA 

C7H8N2O2 

100-15-2 

l-Methyl-2,4,5-trinitroimidazole 

MTNI 

C4H3N5O6 

19183-20-1 

3-Nitro-l,2,4-triazol-5-one 

NTO 

C2H2N4O3 

932-64-9 

1, 3, 5-Triamino-2, 4,6- 
trinitrobenzene 

TATB 

C6H6N6O6 

3058-38-6 

2.4.6.8.10.12- hexanitro- 

2.4.5.8.10.12- 
hexaazaisowurtzitane 

S-CL20 

Yes 

C6H6N12O12 

135285-90-4 

l,3,5-Trinitro-l,3,5- 

triazacyclohexane 

RDX 

Yes 

C3H6N6O6 

121-82-4 

1, 3,5, 7-Tetranitro-l, 3,5,7- 
tetrazacyclooctane 

HMX 

Yes 

C4H8N8O8 

2691-41-0 

N02 


RDX 


HMX 

Figure  1.  Structure  of  compounds  used  for  validation. 
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Figure  2.  Structures  of  compounds  of  interest  in  this  work. 


TURBOMOLE  computations 

The  molecular  structure  was  used  to  generate  a  mapping  of  atoms  in  each 
molecule  in  a  Cartesian  coordinate  system.  This  mapping  used  several 
programs  available  as  freeware/shareware: 

•  Chemdraw  -  for  drawing  molecules  and  generating  SMILES  code 
(Simplified  Molecular  Input  Line  Entry  Specification 

•  Corina  -  for  conversion  of  SMILES  code  to  a  PDB  (protein  databank) 
format  file,  containing  3D  coordinates  of  all  atoms  and  bonding  infor¬ 
mation 

•  Babel  for  conversion  of  the  PDB  format  file  to  XYZ  format  (listing  of 
Cartesian  coordinates  for  each  atom  in  the  compound). 


The  XYZ  format  file  is  then  converted  to  the  required  format  for 
TURBOMOLE  using  a  utility  provided  with  the  TURBOMOLE  software.  A 
sample  of  the  PDB  and  XYZ  file  contents  are  shown  in  Appendix  A  for  the 
reference  compound,  S-CL20. 


For  each  compound,  density  functional  theory  (DFT)  calculations  were 
performed  using  TURBOMOLE  (Version  5.9)  (Ahlrichs  et  al.  1989, 
Trentler  and  Ahlrichs  1995,  Sierka  et  al.  2003).  The  DFT  computation  re¬ 
sulted  in  an  optimized  structure  for  each  compound.  A  second  DFT  com¬ 
putation  was  performed  to  yield  the  input  geometry  and  energy  required 
as  input  to  the  COSMOtherm  program.  The  TZVP  (triple  zeta  valence 
Gaussian)  basis  set  was  employed  (Schafer  et  al.  1994).  These  computa¬ 
tions  were  performed  on  a  Next  Dimension  Next  Demand-II  computer 
system  equipped  with  dual  Advanced  Micro  Devices,  Inc.  (AMD)  Opteron 
processors.  COSMOtherm  (Version  1.06)  was  used  to  predict  physico¬ 
chemical  properties. 
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TATB  MTNI  DNP 

DNAN  MNA  NTO 


Figure  3.  Electrostatic  potential  surfaces  of  compounds  examined  in  this  work. 

The  *. cosmo  file  used  as  input  for  COSMOtherm  is  generated  using  the 
TURBOMOLE  software.  Examples  of  these  input  files  will  be  provided  elec¬ 
tronically  on  request.  The  contents  of  the  *.cosmo  file  provide  a  mapping  of 
the  electrostatic  potential  on  the  surface  of  the  molecule.  Figure  3  shows 
these  surfaces  for  the  reference  compounds  and  the  compounds  of  interest. 

Estimation  of  physical  properties 

The  development  of  estimation  techniques  for  physical  and  thermodynam¬ 
ic  properties  often  relies  on  data  specific  to  the  substance  of  interest.  Often 
used  in  correlation  are  the  critical  point  of  a  substance  (critical  tempera¬ 
ture,  critical  pressure),  the  normal  boiling  point,  the  acentric  factor,  the 
melting  point  and  enthalpies  of  fusion  and  vaporization,  and  the  vapor 
pressure  of  the  pure  component.  For  the  compounds  of  interest  in  this 
work,  very  limited  data  were  available  for  these  physical  properties.  The 
following  paragraphs  review  means  to  estimate  values  for  these  properties 
from  readily  available  information  such  as  molecular  structure. 

Estimation  of  physical  properties  for  the  compounds  of  interest  and  the 
reference  compounds  was  undertaken.  Experimental  data  for  these  physi¬ 
cal  properties  are  sparse.  These  limited  data  were  used,  where  available,  to 
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validate  the  predictive  approach.  The  program,  Cranium  (Molecular 
Knowledge,  Inc.,  Version  1.0.3)  was  used  to  predict  critical  properties 
(temperature,  pressure,  volume  compressibility),  the  acentric  factor,  nor¬ 
mal  boiling  point,  melting  point,  heat  of  formation,  free  energy  of  for¬ 
mation,  heat  of  fusion  at  the  melting  point,  and  latent  heat  of  vaporization 
at  the  normal  boiling  point. 

Critical  properties 

Using  the  molecular  structure  to  estimate  critical  properties  is  the  basis  of 
group  contribution  methods.  These  methods,  which  were  used  to  estimate 
critical  temperature,  critical  pressure  and  critical  volume,  are  based  on 
correlative  approaches  using  experimental  data  coupled  with  regression  to 
obtain  the  contributions  to  a  given  property  from  a  particular  functional 
group.  The  widely  used  Joback  method  (Poling  et  al.  2001),  which  has 
been  the  historical  method  of  choice,  includes  contributions  for  all  func¬ 
tional  groups  of  interest  in  the  present  work.  The  Joback  method  uses  the 
normal  boiling  point  and  contributions  from  the  representative  functional 
groups  to  estimate  critical  temperature.  In  recent  years,  a  number  of  alter¬ 
nate  group  contribution  methods  have  been  proposed  for  estimation  of 
critical  properties;  these  methods  are  based  on  more  extensive  databases, 
in  attempts  to  address  issues  associated  with  first  order  group  contribu¬ 
tions  methods  such  as  failure  to  distinguish  between  isomers.  The  method 
of  Constantinou  and  Gani  (1994)  is  a  second  order  group  contribution 
method  that  allows  isomers  to  be  distinguished  through  the  inclusion  of 
the  second  order  groupings.  However,  in  this  method,  contributions  for 
nitro  groups  (-N02)  bonded  to  a  nitrogen  atom  in  a  ring  are  not  included; 
thus,  this  method  will  not  be  applicable  for  RDX  and  HMX,  which  were 
chosen  as  validation  molecules  in  this  work.  The  method  of  Valderrama 
and  Alvarez  (2006)  was  developed  to  be  applicable  to  high  molecular 
weight  compounds,  such  as  large,  multifunctional  compounds  and  bio¬ 
molecules.  The  Joback  method  and  the  Valderrama/Alvarez  methods  were 
evaluated  in  the  present  effort  for  critical  property  evaluation. 

In  a  previous  effort  (Boddu  et  al.  2004),  the  Joback  method  (Joback  and 
Reid  1987)  was  used  to  estimate  critical  temperature  and  critical  pressure 
for  the  three  reference  compounds,  with  normal  boiling  point  estimated 
with  the  Joback  method  as  modified  by  Stein  and  Brown  (1994).  The 
Joback  method  uses  the  following  expressions  to  evaluate  critical  tempera¬ 
ture,  Tc,  critical  pressure,  Pc,  and  critical  volume,  Vc: 
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(K)  =  Tb 


0.584  +  0.965 


I 


Nk(tck) 


Eq  1 


Pc  (bar)  = 


0.1 13  +  0.0032  *  Natoms 


Nk(pck) 


-2 


Eq  2 


Vc  (cm3  /  mol)  = 


17.5 


+  ZNk(vck) 


Eq  3 


where: 

Nk 

tck,  pck  and  vck 

Natoms 

Tb 


the  number  of  functional  groups  of  type  k  in  the  molecule 
the  contributions  to  the  critical  temperature,  critical  pres¬ 
sure,  and  critical  volume  from  functional  group  k 
the  total  number  of  atoms  in  the  molecule 
the  normal  boiling  point  of  the  substance,  and  can  be  es¬ 
timated  if  an  experimental  value  is  not  available. 


Values  for  the  critical  temperature  and  pressure  were  also  estimated  using 
the  method  of  Valderrama  and  Alvarez  (2006).  This  method  does  not  re¬ 
quire  use  of  the  normal  boiling  point  in  its  estimation  of  critical  tempera¬ 
ture.  It  is  based  on  the  “idea  of  group  contributions  applied  to  estimate  the 
force  and  the  volume  constants  of  the  van  der  Waals  equation  of  state”  and 
“could  be  applied  to  large  complex  substances  as  those  present  in  many 
natural  products”  (Valderrama  and  Alvarez  2006).  Functional  group  contri¬ 
butions  for  the  quantities  Tc  /  ^/p” ,  Tc  /  Pc ,  and  Vc  are  given  by  Valderrama 

and  Alvarez.  These  are  combined  using  the  following  expressions: 


=  dj 


Vc=a3+  XnfVc] 


P3, 


Eq  4 


Eq  5 


Eq  6 


where: 

ni  is  the  number  of  functional  groups  of  type  i  in  the  molecule 
(To/V^  )i’(Tc/Pc  )i>and(V°)  i  are  the  contributions  to  the  respective 

functionals  from  functional  group  i 
ai  and  bi  are  specified  constants. 
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Table  2  lists  the  estimated  critical  temperature  and  pressure  from  both 
methods.  Experimental  critical  temperatures  were  available  for  HMX  and 
RDX.  Using  the  Stein-Brown  modification  of  the  Joback  method,  the  criti¬ 
cal  temperature  for  RDX  was  estimated  at  842.7  K  compared  to  the  exper¬ 
imental  value  of  840  K  (Maksimov  1992);  for  HMX,  the  estimated  critical 
temperature  was  913.5  K  compared  to  the  experimental  value  of  927  K 
(Maksimov  1992).  The  method  of  Valderrama  and  Alvarez  estimated  criti¬ 
cal  temperatures  for  RDX  and  HMX  that  were  much  greater  than  the  val¬ 
ues  reported  in  the  literature:  for  RDX,  1112.3  K;  and  for  HMX  1234.4  K. 


The  values  estimated  using  the  Joback  method  as  modified  by  Stein  and 
Brown  were  in  better  agreement  with  the  available  literature  data  for  the 
reference  compounds.  The  values  for  critical  temperature  and  pressure  es¬ 
timated  using  the  Stein-Brown  modification  of  the  Joback  method  were 
used  in  further  computations. 

Normal  boiling  point  and  acentric  factor 

The  estimation  of  the  normal  boiling  point  and  the  acentric  factor  are  also 
of  interest  in  this  work.  The  Stein-Brown  method  (1994)  was  used  to  esti¬ 
mate  the  normal  boiling  point  of  both  reference  compounds  and  com¬ 
pounds  of  interest  in  this  work.  This  method  uses  the  same  functional 
groups  as  the  Joback  method,  but  applies  a  correction  to  estimated  boiling 
temperatures  above  800  K.  This  correction  is  in  essence  responsible  for 
the  difference  observed  in  the  predicted  critical  temperatures  listed  in  Ta¬ 
ble  2.  The  acentric  factor  is  estimated  using  its  definition,  along  with  the 
estimated  critical  temperature,  critical  pressure  and  normal  boiling  point. 
Table  3  lists  predictions  of  normal  boiling  point  and  acentric  factor  for  the 
compounds  of  interest  and  reference  compounds. 


Table  2.  Critical  properties. 


Literature  Data 

Joback/Stein  Method 

Valderrama/Alveraz  Method 

Compound 

To  (K) 

Tc(K) 

Pc(bar) 

Tc(K) 

Pc(bar) 

DNAN 

806.3 

39.9 

938.9 

34.67 

DNP 

816.8 

69.33 

1036.4 

59.72 

MNA 

747.5 

41.68 

845.8 

41.34 

MTNI 

845.0 

54.71 

1105.8 

37.21 

NTO 

828.5 

86.05 

1291.6 

145.15 

TATB 

912.7 

56.45 

1189.0 

31.14 

S-CL20 

1058.0 

48.90 

1539.4 

21.96 

RDX 

840a 

842.7 

58.01 

1112.3 

35.31 

HMX 

927a 

913.5 

52.97 

1234.4 

28.83 

a  Maksimov  (1992). 
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Table  3.  Estimated  normal  boiling  point  and  acentric  factor. 


Compound 

Estimated  Normal  Boiling  Point  Tnbp  (K) 
Stein-Brown  Method 

Estimated  Acentric  Factor  co 
Definition  of  Acentric  Factor 

DNAN 

587.6 

0.858 

DNP 

574.5 

0.307 

MNA 

527.3 

0.659 

MTNI 

629.1 

0.707 

NTO 

568.4 

0.798 

TATB 

711.6 

1.802 

S-CL20 

861.8 

2.547 

RDX 

626.6 

1.252 

HMX 

709.5 

1.700 

Vapor  pressure 

There  are  very  few  experimental  data  available  for  the  vapor  pressure 
(over  liquid  or  solid  phases)  for  the  compounds  of  interest  and  the  refer¬ 
ence  compounds.  Rosen  and  Dickenson  (1969)  reported  on  vapor  pressure 
and  heat  of  sublimation  for  high  melting  organic  explosives,  including  ref¬ 
erence  compounds  RDX  and  HMX,  and  one  compound  of  interest,  TATB. 

Taylor  and  Crookes  (1976)  also  reported  on  these  properties  for  HMX. 
Cundall  et  al.  (1978)  provided  measurements  of  vapor  pressure  for  a  slate 
of  compounds  including  RDX  and  HMX.  Dionne  et  al.  (1986)  measured 
the  vapor  pressure  of  RDX.  Sinditskii  et  al.  (2006)  reported  the  vapor 
pressure  over  both  liquid  and  solid  phase  NTO,  fitting  their  data  to  a  sim¬ 
ple  Claperyon  relationship.  These  data  were  used,  where  possible,  for 
comparison  with  predictions  of  vapor  pressure  for  the  compounds  of  in¬ 
terest  and  the  reference  compounds. 

COSMOtherm  was  used  to  predict  the  vapor  pressure  for  each  reference 
compound  and  for  the  compounds  of  interest.  Table  4  lists  Antoine  coeffi¬ 
cients  for  the  COSMOtherm  vapor  pressure  predictions.  For  those  com¬ 
pounds  where  experimental  data  were  available  in  the  literature,  a  com¬ 
parison  was  made  with  COSMOtherm  predictions.  For  RDX,  literature 
data  were  available  from  three  independent  sources.  Figure  4  shows  these 
data  with  the  COSMOtherm  prediction.  The  COSMOtherm  predictions  are 
much  larger  than  the  reported  experimental  vapor  pressures,  and  there  is 
considerable  scatter  among  the  three  literature  data  sets.  The  scatter 
among  these  measurements  is  indicative  of  the  difficulty  in  measuring  va¬ 
por  pressure  for  compounds  exhibiting  very  low  volatility. 
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Table  4.  Antoine  coefficients  generated  using  COSMOtherm. 


Compound 

Aj 

Bj 

Ci 

DNAN 

18.01 

7087.5 

-43.94 

DNP 

17.72 

7092.4 

-60.25 

MNA 

18.19 

6757.9 

-61.43 

MTNI 

17.99 

7080.5 

-41.10 

NTO 

20.63 

7334.04 

-76.29 

TATB 

18.27 

6756.2 

-62.45 

S-CL20 

18.74 

9247.9 

-38.79 

RDX 

18.56 

7718.5 

-42.65 

HMX 

20.24 

11881.4 

-35.96 

Where  ln[P(mbar)]=Ai-[Bi/[T(K)+Ci]] 

Reciprocal  Temperature,  in,  [1JK] 

Figure  4.  Vapor  pressure  of  RDX  shown  with  literature  data. 

Figure  5  shows  HMX  with  experimental  data  available  in  three  independ¬ 
ent  sources.  Again,  COSMOtherm  predicts  much  larger  vapor  pressure  at  a 
given  temperature  compared  to  the  experimental  data.  Figure  6  shows  the 
vapor  pressure  of  TATB  as  measured  by  Rosen  and  Dickenson  (1969) 
against  the  COSMOtherm  prediction. 


Vapor  Pressure,  Psat,  (Pa)  vapor  Pressurei  Psat,  pa) 
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Reciprocal  Temperature,  1/T,  (1/K) 

Figure  5.  Vapor  pressure  of  HMX  shown  with  literature  data. 


Reciprocal  Temperature,  1/T,  (1/K) 

Figure  6.  Vapor  pressure  of  TATB  shown  with  literature  data. 
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Figure  7.  COSMOtherm  predictions  of  vapor  pressure  for  all  compounds 
examined  in  this  work. 

Again,  the  experimental  data  from  the  literature  are  much  smaller  than  the 
COSMOtherm  predictions.  Figure  7  shows  COSMOtherm  predictions  for 
all  compounds  of  interest  and  for  the  reference  compounds.  COSMOtherm 
predicts  the  trend  in  decreasing  volatility  to  be: 

MTNI  >  DNAN  >  MNA  >  TATB  >  RDX  >  DNP  >  NTO  »  e-CL20  »  HMX. 

Melting  point  and  enthalpy  of  fusion 

The  melting  point  for  each  of  the  compounds  of  interest  and  the  reference 
compounds  were  available  from  the  literature.  Literature  enthalpy  of  fu¬ 
sion  data  were  available  for  only  three  compounds  (RDX,  NTO,  and 
s— CL20).  Domalski  and  Hearing  (1996)  reported  the  enthalpy  of  fusion  of 
RDX  as  37.66  kJ/gmol.  For  NTO,  the  reported  values  of  27  kJ/mol  (Liu  et 
al.  1995)  and  92.8  kJ/mol  (Kim  et  al.  1998)  differ  by  a  factor  of  4.  For 
8-CL20,  the  reported  enthalpy  of  fusion  value  was  42.7  kJ/mol,  estimated 
from  data  for  similar  nitro  compounds  (Zeman  and  Jalovy  2000);  it  was 
not  an  experimentally  determined  value. 

The  method  of  Chickos  et  al.  (1999)  was  employed  for  estimation  of  the 
total  phase  change  entropy  at  the  melting  point.  In  this  group  contribution 
method,  the  total  entropy  of  fusion  becomes  numerically  equivalent  to  the 


x  MTNI 
-DNAN 
o  MNA 
&  TATB 

□  RDX 

+  DNP-34 
a  NTO 

□  CL20 
o  HMX 
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entropy  of  fusion  when  there  is  a  single  solid-liquid  phase  transition.  En¬ 
thalpy  of  fusion  was  then  computed  from  the  predicted  entropy  of  fusion 
and  the  literature  melting  point  for  each  compound.  Table  5  summarizes 
the  estimated  enthalpy  of  fusion  for  the  compounds  in  this  work,  and  (if 
available)  the  literature  data  for  melting  point  and  enthalpy  of  fusion. 

The  estimated  enthalpy  of  fusion  for  RDX  at  39.52  kJ/mol  is  in  good 
agreement  with  the  experimental  value  of  37.66  (Domalski  and  Hearing 
1996).  For  8-CL20,  the  estimated  value  is  much  lower  at  13.7  kJ/mol  com¬ 
pared  to  the  literature  value  of  42.7  (Zeman  and  Jalovy  2000).  For  NTO, 
the  estimated  value  of  20.17  kJ/mol  was  comparable  to  the  reported  litera¬ 
ture  value  of  27  kJ/mol.  The  estimated  enthalpy  of  fusion  values  were  used 
in  all  further  predictive  efforts. 


Table  5.  Melting  point  and  enthalpy  of  fusion. 


Compound 

Melting  Point 

Tm  (K) 
(reference) 

Estimated  Enthalpy 
of  Fusion, 

AHm  (kJ/mol) 

Literature  Enthalpy 
of  Fusion, 

AHm  (kJ/mol) 
(reference) 

DNAN 

359. 9a 

19.91 

DNP 

363-364b 

27.6 

MNA 

425.15° 

18.85 

MTNI 

355. 15d 

19.82 

NTO 

539. 35e,  547.9f 

20.17 

27e,  92. 8f 

TATB 

594-599^ 

43.0 

S-CL20 

513h 

13.7 

42.7J 

RDX 

478.5' 

39.52 

37.66k 

HMX 

553. 15a 

45.0 

a  Maksimov  (1992). 

b  Beringer  et  al.  (1953). 

c  Katritzky  et  al.  (2005),  p  179. 

d  Material  Safety  Data  Sheet  [OSHA]  (MSDS)  for  N-methyl-4-nitroaniline, 

htto://otcl. chem.ox.ac.uk/MSDS/ME/N-methvl-4-nitroaline.html .  accessed 

01/05/07. 

e  Liu  et  al.  (1995). 

f  Kim  et  al.  (1998). 

gChoetal.  (2002). 

h  Andelkovic  (2000). 

1  Hall  (1971). 

j  Zeman  and  Jalovy  (2000). 

k  Domalski  and  Hearing  (1996). 
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Table  6.  Abraham  descriptors  for  compounds  of  interest  and  reference  compounds. 


Compound 

Excess  Molar 
Refraction 

r2 

Polarity/ 

Polarizability 

71 2 

Hydrogen 
Bond  Acidity 

Ia2 

Hydrogen 

Bond 

Basicity 

Ip? 

McGowan 

Characteristic 

Volume 

Vx 

DNAN 

1.7059 

1.9993 

0.0057 

0.6140 

1.3289 

DNP 

1.7891 

1.8702 

1.2048 

0.3796 

0.9655 

MNA 

0.9795 

1.4126 

1.4126 

0.6339 

1.1521 

MTNI 

3.5180 

3.0367 

0.0845 

0.4823 

1.2902 

NTO 

0.9991 

1.5359 

1.5333 

0.5891 

0.8093 

TATB 

1.8381 

2.3187 

0.3112 

0.8470 

1.4534 

S-CL20 

7.0001 

5.1712 

1.0047 

0.3969 

2.1330 

RDX 

3.0962 

2.9028 

0.2393 

0.5430 

1.2475 

HMX 

4.0232 

3.8679 

0.6605 

0.7576 

1.5855 

Computation  of  molecular  descriptors 

In  addition  to  the  conventional  physical  properties  estimated  from  molec¬ 
ular  structure,  the  Abraham  descriptors  were  computed.  These  descriptors 
are  the  foundation  of  the  Abraham  method  (Abraham  1993)  that  correlates 
solvent-solute  interactions  using  a  solvation  parameter  approach.  The 
solvation  property  (SP)  is  a  linear  combination  of  solute  specific  de¬ 
scriptors,  as  shown  in  equation  (7): 

log(SP)  =  c  +  rR2  +S7I2  +a^a“+b^p2  +vVx  ? 

where: 

=  the  excess  molar  refraction  (molar  refraction  of  the  solute  minus 
the  molar  refraction  of  an  alkane  of  equivalent  volume 
=  the  descriptor  for  combined  dipolarity/polarizability 
=  the  overall  solute  hydrogen  bond  acidity 

=  the  overall  solute  hydrogen  bond  basicity 

=  McGowan’s  characteristic  volume  (Abraham  and  McGowan  1987). 

Table  6  lists  the  Abraham  descriptors  as  computed  in  COSMOtherm. 

Estimation  of  physicochemical  properties 

Physicochemical  properties  include  aqueous  solubility,  octanol-water  par¬ 
tition  coefficient,  and  Henry’s  constant  (among  others).  These  properties 
determine  how  a  particular  compound  partitions  in  the  environment  and 
distributes  among  phases.  Key  to  these  properties  is  an  understanding  of 
their  influence  on  potential  treatment  technologies.  For  example,  alkaline 


R2 


I«? 

V. 
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hydrolysis  is  often  used  to  remediate  contaminated  water  streams,  but  the 
solubility  of  the  contaminant  in  the  stream  is  critical;  alkaline  hydrolysis 
may  not  be  the  preferred  method  of  remediation  for  compounds  exhibiting 
extremely  low  solubilities  in  aqueous  solution. 

Aqueous  solubility 

Monteil-Rivera  et  al.  (2004)  reported  experimental  measurements  for  the 
aqueous  solubility  as  a  function  of  temperature,  the  i-octanol/water  parti¬ 
tion  coefficient  at  21  °C,  and  the  soil-water  distribution  coefficient  at  21  °C, 
for  s  -CL20,  RDX  and  HMX.  Aqueous  solubility  for  s  -CL20  has  also  been 
reported  by  Karakaya  et  al.  (2005).  Lynch  et  al.  (2001)  measured  the  solu¬ 
bility  of  three  energetic  compounds,  trinitrotoluene  (TNT),  RDX  and  HMX 
over  a  moderate  temperature  range  (2-38  °C)  as  a  function  of  solution  pH 
(4.2-6.2).  Aqueous  solubility  increased  with  temperature  increase.  Monteil- 
Rivera  et  al.  (2004)  also  reported  aqueous  solubility  measurements  for 
RDX,  HMX,  and  S-CL20  over  the  temperature  range  of  5  to  60  °C.  These 
data  provide  the  means  for  validating  the  estimation  techniques. 

For  the  compounds  of  interest,  experimental  aqueous  solubility  data  have 
only  been  reported  for  NTO,  MNA,  and  DNAN.  Kim  et  al.  (1998)  reported 
measurements  of  aqueous  solubility  of  NTO  in  water  over  a  fairly  wide 
temperature  range  of  284.65  to  367.55  K.  Boddu  et  al.  (2008a)  reported 
experimental  measurements  of  aqueous  solubility  for  DNAN.  Experi¬ 
mental  solubility  measurements  for  MNA  were  also  reported  by  Boddu  et 
al.  (2008b). 


Group  contribution  methods  for  the  estimation  of  aqueous  solubility  at 
25  °C  are  numerous  and  vary  in  their  complexity.  The  methods  applicable 
in  this  work  are  limited  as  many  of  the  reported  methods  do  not  include  a 
functional  group  for  nitrogen  in  a  ring  structure  in  their  development.  The 
widely-used  methods  of  Klopman  (1992)  contain  this  required  functional 
group  in  Model  II  of  Klopman.  In  Klopman’s  Model  II,  67  different  func¬ 
tional  groups  are  considered.  The  general  form  of  the  Klopman  model  is: 


log 


10 


=  c,  +  £gA 


Eq  8 


where: 


Sw 

Co 

Gi 

gi 


the  aqueous  solubility  in  g/L 
a  constant 

the  contribution  to  the  solubility  due  to  functional  group  i 
the  number  of  functional  group  i  contained  in  the  molecule. 
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The  solubility  of  a  solid  solute  in  a  liquid  is  given  by  (Poling  et  al.  2001): 


ln(y2x2)  =  - 


AS„ 


R 


(  t  'N 

L  m  _  j 

It  J 


AH. 


RT 


m  J 


Eq  9 


where: 

X2  =  the  mole  fraction  of  the  solute  in  the  solution 
ASm  =  the  entropy  of  fusion,  AHm  is  the  enthalpy  of  fusion 

Tm  =  the  melting  point  in  K 
R  =  the  gas  constant 
T  =  the  temperature  of  interest  in  K. 

The  ideal  solubility  is  obtained  when  y2,  the  activity  coefficient  of  the  so¬ 
lute  in  the  solution,  is  set  to  unity. 

Table  7  lists  a  compilation  of  aqueous  solubility  predictions  using  the 
method  of  Klopman  et  al.,  COSMOtherm,  and  the  ideal  solubility  (equa¬ 
tion  [8]).  The  predictions  using  the  method  of  Klopman  and  from 
COSMOtherm  are  in  excellent  agreement  with  one  another  for  every  com¬ 
pound  examined.  The  ideal  solubility  predicted  for  each  compound  is 
much  greater  (often  3  to  4  orders  of  magnitude  greater)  compared  to  pre¬ 
dictions  of  Klopman  and  of  COSMOtherm.  These  results  indicate  that 
aqueous  solutions  containing  these  energetic  compounds  exhibit  extreme¬ 
ly  large  deviations  from  ideal  behavior. 

Also  of  interest  in  the  present  work  is  the  temperature  dependence  of 
aqueous  solubility.  No  group  contribution  methods  exist  to  predict  said 
temperature  dependence.  Thus,  an  alternate  approach  was  undertaken  us¬ 
ing  COSMOtherm.  Aqueous  solubility  data  as  a  function  of  temperature 
were  available  for  all  three  reference  compounds.  Predictions  for  these 
provide  a  means  to  assess  estimation  of  temperature  dependent  aqueous 
solubility  and  were  compared  with  available  literature  data  where  possible. 

Aqueous  solubility  over  the  temperature  range  of  273.15  K  to  373.15  K  was 
estimated  using  COSMOtherm  for  each  of  the  compounds  of  interest  and 
for  the  reference  compounds.  For  both  RDX  and  HMX,  the  aqueous  solu¬ 
bilities  predicted  using  COSMOtherm  were  smaller  than  the  literature  data 
(Figures  8  and  9,  RDX  and  HMX,  respectively).  The  experimental  data  of 
Lynch  et  al.  (2001)  and  of  Monteil-Rivera  et  al.  (2004)  are  in  excellent 
agreement.  Figure  10  shows  aqueous  solubility  data  for  8-CL20;  the  litera¬ 
ture  data  of  Monteil-Rivera  et  al.  (2004)  and  of  Karakaya  et  al.  (2005)  are 
shown  as  are  the  COSMOtherm  predictions.  The  data  of  Karakaya  et  al. 
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(2005)  extend  over  a  fairly  broad  temperature  range  (273  -  338  K).  The 
predictions  from  COSMOtherm  reflect  the  general  trends  exhibited  by  the 
experimental  data. 


Table  7.  Comparison  of  predictive  methods  for  aqueous  solubility  at  298.15  0  K. 


Compound 

Aqueous  Solubility,  S 
(g/L)  from  literature 

Aqueous  Solubility,  S 
(g/L) 

Method  of  Klopman 
etal.  (1992) 

Aqueous  Solubility,  S 
(g/L) 

COSMOtherm 

prediction 

Ideal  Solubility,  S 
(g/L) 

Equation  (8) 

DNAN 

0.2769 

0.14 

0.159 

3706.5 

DNP 

0.51 

0.733 

1390.7 

MNA 

0.085b 

1.17 

0.347 

971.3 

MTNI 

0.23 

0.291 

4619.9 

NTO 

9.99  at  11.5  °Cc 

19.99  at  33.9  °Cc 

7.56 

9.97 

55.8 

TATB 

0.16 

0.197 

2.5 

S-CL20 

3.65  x  10-3  d 

3.86  x  lO  3 

2.91  x  lO  3 

17.9 

RDX 

56.35  x  lO'3  d 

14.2  x  lO'3 

12.6  x  lO  3 

40.3 

HMX 

4.46  x  IQ  3  d 

2.7  x  lO  3 

2.16  x  lO  3 

3.81 

a  Boddu  et  al.  (2008a). 

bBodduetal.  (2008b). 

c  Kim  etal.  (1998). 

d  Monteil-Rivera  et  al.  (2004). 
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Figure  8.  Aqueous  solubility  of  RDX;  literature  data  shown  with 
COSMOtherm  predictions. 
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Figure  9.  Aqueous  solubility  of  HMX;  literature  data 
shown  with  COSMOtherm  prediction. 
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Figure  10.  Aqueous  solubility  of  e-CL20;  literature  data 
shown  with  COSMOtherm  predictions. 
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Figure  n  shows  the  experimental  data  of  Kim  et  al.  (1998)  for  the  aqueous 
solubility  of  NTO.  Again,  the  general  trend  of  the  experimental  data  is  well 
predicted  by  COSMOtherm.  The  COSMOtherm  predictions  are  slightly 
lower  than  the  experimental  data. 

Figures  12  and  13  provide  a  comparison  of  recently  reported  experimental 
data  for  DNAN  (Boddu  et  al.  2008a)  and  MNA  (Boddu  et  al.  2008b),  re¬ 
spectively.  The  COSMOtherm  predictions  for  DNAN  reflect  the  appropri¬ 
ate  temperature  dependence  over  the  small  temperature  range  examined 
(298.15  -  318.15  K)  compared  to  the  literature  data,  but  predictions  are 
slightly  lower  than  the  experimental  data.  For  MNA,  the  experimental 
aqueous  solubility  is  smaller,  and  COSMOtherm  predictions  are  in  good 
agreement  with  the  experimental  data.  At  298.15  K,  the  experimental  sol¬ 
ubility  for  DNAN  is  0.276  g/L  while  COSMOtherm  predicts  0.159  g/L.  At 
298.15  K,  the  experimental  solubility  for  MNA  is  0.085  g/L  while  the 
COSMOtherm  prediction  is  slightly  larger  at  0.347  g/L. 

In  general,  COSMOtherm  predictions  of  aqueous  solubility  are  quite  good, 
considering  that  these  are  a  priori  estimates  based  on  molecular  structure 
and  on  physical  properties  predicted  using  a  variety  of  group  contribution 
methods.  COSMOtherm  predictions  were  within  10  percent  of  the  experi¬ 
mental  data  values,  where  available. 

Figure  14  shows  COSMOtherm  predictions  for  all  compounds  of  interest 
and  for  the  reference  compounds.  The  aqueous  solubility  for  all  com¬ 
pounds  shows  strong  temperature  dependence.  NTO  is  predicted  to  be  ex¬ 
tremely  soluble  in  water,  with  solubilities  on  the  order  of  ~  10  g/L  at  ambi¬ 
ent  temperature.  The  other  compounds  of  interest  are  all  predicted  to  be 
more  soluble  than  the  three  reference  compounds,  by  as  much  as  3  orders 
of  magnitude. 

Octanol-water  partition  coefficient 

The  octanol-water  partition  coefficient  provides  a  measure  of  the  hydro- 
phobicity  of  a  compound  in  addition  to  serving  as  a  measure  of  the  affinity 
of  a  solute  to  solubilize  in  an  organic-rich  (l-octanol)  phase  compared  to  in 
an  aqueous  phase  (water).  The  partition  coefficient  is  often  reported  as 
log  io(Kow),  with  a  positive  value  indicating  that  a  solute  preferentially  par¬ 
titions  into  the  organic  phase,  while  a  negative  value  indicates  that  the  so¬ 
lute  preferentially  partitions  into  the  aqueous  phase.  A  positive  octanol- 
water  partition  coefficient  also  indicates  that  the  particular  solute  would 
have  an  affinity  for  organic  matter  in  soil. 
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Figure  11.  Aqueous  solubility  of  NTO;  literature  data 
shown  with  COSMOtherm  predictions. 


Figure  12.  Aqueous  solubility  of  DNAN;  literature  data 
shown  with  COSMOtherm  predictions. 
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Figure  13.  Aqueous  solubility  of  MNA;  literature 
data  shown  with  COSMOtherm  predictions. 
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Figure  14.  Solubilities  of  compounds  of  interest  and  reference  compounds. 
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Two  group  contribution  methods  were  examined  for  prediction  of  the 
octanol-water  partition  coefficient:  the  method  of  Broto  et  al.  (1984)  and 
the  method  of  Ghose  et  al.  (1988).  These  methods  are  atom  contribution 
methods,  wherein  contributions  from  each  atom  are  combined  with  con¬ 
tributions  dependent  on  hybridization  and  the  presence  of  heteroatoms  as 
neighbors.  These  methods  both  used  extensive  training  sets,  composed  of 
both  a  diversity  of  atoms  and  atomic  environments.  The  Broto  method  is 
reported  to  provide  Kow  values  precise  to  0.4  log  units,  while  the  precision 
of  the  Ghose  method  is  reported  as  0.496  log  units  (Ghose  et  al.  1988). 

COSMOtherm  predictions  for  the  octanol-water  partition  coefficient  were 
also  performed.  Table  8  lists  the  available  experimental  data,  the  predic¬ 
tions  using  the  group  contribution  methods,  and  the  COSMOtherm  predic¬ 
tions  at  298.15  K.  For  the  majority  of  compounds,  the  predictions  from  the 
group  contribution  methods  are  in  agreement.  For  DNP,  the  methods  pro¬ 
vide  estimates  that  are  similar  in  magnitude,  but  opposite  in  sign.  For 
RDX,  the  values  predicted  using  the  group  contribution  methods  and  us¬ 
ing  COSMOtherm  are  in  good  agreement  with  the  available  literature  val¬ 
ues.  The  predictions  for  MNA  using  either  group  contribution  method  or 
COSMOtherm  are  significantly  smaller  than  the  experimentally  deter¬ 
mined  values  of  2.1028  +  0.02  and  2.04.  A  difference  of  1  unit  in 
log  10  (Kow)  translates  into  an  order  of  magnitude  difference  in  the  parti¬ 
tion  coefficient. 

Based  on  the  experimental  data  available  in  the  literature,  DNAN,  MNA, 
and  8-CL20  would  preferentially  partition  in  the  organic  (l-octanol)  phase 
to  a  greater  extent  than  either  RDX  or  HMX.  MTNI  is  predicted  to  also 
partition  into  the  organic  phase.  Both  TATB  and  NTO  are  predicted  to 
preferentially  partition  into  the  aqueous  phase. 

Henry’s  constant 

The  partitioning  of  a  solute  between  the  gas  phase  and  the  aqueous  phase 
is  described  using  the  Henry’s  constant.  This  physicochemical  parameter 
was  evaluated  using  COSMOtherm  results  for  the  aqueous  solubility  and 
the  vapor  pressure  of  the  solute  over  the  temperature  range  273.15  to 
373.15  K.  Figure  15  shows  the  resulting  predictions  for  all  compounds  ex¬ 
amined.  With  the  exception  of  TATB,  the  Henry’s  constant  values  indicate 
that  the  compounds  prefer  to  partition  into  the  aqueous  phase,  rather  than 
into  a  gas  phase. 
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Table  8.  Octanol-water  partition  coefficient  at  298.15  °  K. 


Compound 

Experimental 

Data 

Method  of  Broto 
etal. 

Method  of  Ghose 
etal. 

COSMOtherm 

Prediction 

DNAN 

1.6128 

1.38 

1.70 

1.92 

DNP 

0.9 

-0.99 

0.37 

MNA 

2.1028  +  0.02b 

2.04c 

1.47 

1.57 

0.81 

MTNI 

- 

0.48 

1.64 

NTO 

-0.82 

- 

-1.19 

TATB 

-1.5 

-0.44 

4.74 

8-CL20 

1.92d  at  21  °C 

- 

1.38 

4.14 

RDX 

0.87e 

0.90d  at  21  °C 

0.85 

1.04 

0.90 

HMX 

0.165d  at  21  °C 

1.13 

1.39 

0.42 

a  Boddu  etal.  (2008). 

bBodduetal.  (2008). 

c  Banerjee  et  al.  (1980). 

d  Monteil-Rivera  et  al.  (2004).  Measurements  at  21  °C. 

e  Hansch  (1979). 

Temperature,  T,  (K) 

Figure  15.  Predicted  temperature  dependence  of  Henry’s  Constant. 


ERDC/CERL  TR-10-27 


25 


The  extremely  low  volatility  of  the  compounds  is  the  dominant  factor  in 
determining  Henry’s  constant.  Table  9  lists  the  predictions  for  Henry’s 
constant  using  COSMOtherm  along  with  available  experimental  data  for 
MNA  and  DNAN  (Boddu  et  al.  2008a,  2008b).  The  estimated  values  are 
greater  than  the  predicted  values,  by  at  least  one  order  of  magnitude.  The 
influence  of  vapor  pressure  (estimated  using  COSMOtherm)  on  the  pre¬ 
dicted  values  warrants  investigation.  Accurate  measurement  of  vapor 
pressure  for  the  compounds  of  interest  over  the  temperature  range  283  K 
to  333  K  would  provide  for  improved  estimates  of  the  Henry’s  constant. 


Table  9.  Henry’s  constant  at  298.15  0  K. 


Compound 

Experimental  Data 
(m3-Pa/mol) 

COSMOtherm 

Prediction 

(m3-Pa/mol) 

DNAN 

1.3669 

6.37  x  lO-3 

DNP 

1.20  x  lO  5 

MNA 

0.61b 

6.18  x  lO  3 

MTNI 

5.26  x  lO  3 

NTO 

4.45  x  lO  7 

TATB 

0.397 

S-CL20 

6.73  x  lO  4 

RDX 

0.153 

HMX 

1.75  x  lO  7 

a  Boddu  et  al.  (2008). 

b  Boddu  etal.  (2008). 
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3  Discussion 

This  work  completed  the  estimation  of  physical  properties  and  physico¬ 
chemical  properties  for  the  compounds  of  interest  and  for  the  reference 
compounds.  One  of  the  critical  variables  in  the  predictive  approach  is  the 
vapor  pressure  (or  sublimation  pressure).  Predictions  of  vapor  pressure  by 
COSMOtherm  were  compared  with  experimental  data  where  available.  For 
RDX,  HMX,  and  TATB,  the  estimated  vapor  pressures  were  larger  than  the 
literature  data  by  at  least  two  orders  of  magnitude.  Even  where  multiple 
data  sets  are  available  (for  example,  RDX,  Figure  4;  HMX,  Figure  5),  there 
are  significant  discrepancies  between  the  data  sets  that  cannot  be  attribut¬ 
ed  to  experimental  error.  The  influence  of  vapor/sublimation  pressure  on 
the  physicochemical  properties,  in  particular,  on  Henry’s  constant  is  sig¬ 
nificant.  One  recommendation  to  improve  the  quality  of  predictions  is  to 
expand  the  existing  experimental  program  at  CERL  to  measure  the  va¬ 
por/sublimation  pressure  of  the  energetic  compounds  of  interest  and  the 
reference  compounds.  Establishing  an  experimental  protocol  for  the  accu¬ 
rate  measurement  of  the  vapor/sublimation  pressures  of  materials  with 
extremely  low  volatility  would  be  of  great  benefit  to  the  present  work. 

Results  for  the  aqueous  solubility  as  a  function  of  temperature  for  the  ref¬ 
erence  compounds  provide  a  measure  of  validation  for  the  predictive  ap¬ 
proach  undertaken.  With  this  approach,  it  is  possible  to  predict  aqueous 
solubility  to  within  10  percent  of  the  experimental  data  using  the  molecu¬ 
lar  structure  alone.  The  influence  of  the  vapor/sublimation  pressure  on 
aqueous  solubility  should  be  examined  as  more  accurate  data  become 
available  for  this  property.  An  examination  of  the  influence  on  aqueous 
solubility  of  salt  addition,  or  the  presence  of  caustic  in  the  aqueous  solu¬ 
tion  would  also  be  beneficial,  as  data  of  this  nature  could  be  used  to  fur¬ 
ther  enhance  predictive  capabilities  using  molecular  structure  alone. 

The  efforts  at  CERL  undertaken  in  recent  years  to  measure  the  physico¬ 
chemical  properties  as  a  function  of  temperature  for  these  energetic  com¬ 
pounds  of  interest  will  provide  additional  experimental  data  that  may  be 
used  to  validate  and  further  refine  the  predictive  techniques.  Boddu  et  al. 
(2008a,  2008b)  reported  the  influence  of  various  inorganic  salts  on  the 
solubility  of  DNAN  and  MNA;  a  demonstration  of  predictive  capabilities 
that  can  reproduce  the  effects  of  salt  addition  on  aqueous  solubility  would 
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be  beneficial.  Such  predictive  capabilities  could  then  be  used  to  exam¬ 
ine/estimate  the  behavior  of  energetic  compounds  developed  in  the  future. 
Data  for  the  partition  coefficients  may  also  be  used  to  examine  whether 
other  techniques  provide  more  accurate  predictions.  It  is  difficult  to  assess 
the  accuracy  of  the  current  predictive  approach  with  the  sparse  amount  of 
experimental  data  available  at  the  present  time,  in  particular  for  Henry’s 
constant. 
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4  Treatment  Technologies 

A  variety  of  treatment  technologies  have  been  examined  for  the  remediation 
of  waste  streams  containing  energetic  compounds.  Rodgers  and  Bunce 
(2001)  published  a  review  of  the  major  technologies  employed  for  remedia¬ 
tion  of  nitroaromatic  explosives.  The  majority  of  published  studies  have  ex¬ 
amined  the  remediation  of  waste  streams  containing  TNT  or  TNT-related 
compounds.  A  few  have  examined  remediation  of  RDX-  and  HMX- 
contaminated  streams,  and  still  fewer  have  examined  treatment  of  streams 
containing  S-CL20.  The  primary  focus  of  this  effort  is  on  waste  water  streams 
containing  energetic  compounds,  which  are  generated  during  production  of 
the  compounds  of  interest.  These  streams  may  contain  an  energetic  com¬ 
pound  at  or  below  its  solubility  limit,  and  may  also  contain  small  amounts  of 
the  compound  in  the  solid  phase.  In  the  present  effort,  it  is  assumed  that  only 
a  single  contaminant  is  present  in  a  given  waste  water  stream. 

A  stream  may  be  remediated  by  the  destruction  of  the  contaminant  into 
less  problematic  compounds;  technologies  resulting  in  destruction  of  the 
contaminant  include  advanced  oxidation  processes,  chemical  reduction 
processes,  and  bioremediation  processes.  Alkaline  hydrolysis  has  been  ex¬ 
amined  as  a  treatment  method  for  streams  contaminated  with  energetic 
materials,  including  RDX,  HMX,  and  8-CL20.  Remediation  may  also  be 
accomplished  with  the  recovery  of  the  contaminant  to  a  sufficient  degree 
from  the  waste  stream  so  that  the  wastewater  may  safely  be  released  to  the 
environment.  The  recovery  is  accomplished  through  separation  processes 
such  as  extraction,  stripping,  adsorption  or  filtration.  In  some  instances,  it 
may  be  beneficial  to  combine  distinct  recovery  and  destruction  steps  in  a 
process  so  as  to  improve  the  overall  efficiency  of  the  remediation  process. 

Hwang  et  al.  (2005)  examined  remediation  of  TNT-contaminated  water  by 
alkaline  hydrolysis.  Batch  reactor  studies  conducted  over  a  range  of  pH 
(11-12)  and  over  a  range  of  initial  contaminant  concentrations 
(5-25  mg/L)  yielded  kinetic  rate  information.  Using  a  pseudo  first-order 
rate  constant  to  describe  the  dependence  on  TNT  concentration,  Hwang  et 
al.  concluded  that  the  destruction  proceeded  at  a  greater  rate  when  the  pH 
was  greater  (i.e.,  more  alkaline).  Balakrishnan  et  al.  (2003)  examined  the 
alkaline  hydrolysis  of  RDX,  HMX,  and  8-CL20  in  two  solutions  of  differing 
pH  (10  versus  12.3).  They  found  that  8-CL20  was  degraded  much  more 
quickly  in  the  alkaline  environment  compared  to  RDX  or  HMX  and  con- 
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eluded  that  this  was  due  to  the  susceptibility  of  the  polycyclic  nitramine 
cage  in  8-CL20  to  nucleophilic  attack. 

Heilman  et  al.  (1996)  also  examined  the  alkaline  hydrolysis  of  RDX  and 
HMX  in  aqueous  solution  and  adsorbed  to  activated  carbon.  They  demon¬ 
strated  that  elevated  temperature  (80  °C)  and  high  pH  (12)  provided  for 
increased  hydrolysis  reaction  rate.  Karakaya  et  al.  (2005)  examined  the 
alkaline  hydrolysis  of  S-CL20  in  both  homogeneous  and  heterogeneous 
environments.  Other  process  variables  examined  included  caustic  concen¬ 
tration  (i.e.,  pH),  and  temperature  (15  to  40  °C).  Karakaya  et  al.  also  found 
that  increased  temperature  and  higher  caustic  loading  (high  pH)  increased 
the  rate  of  degradation  by  alkaline  hydrolysis. 

Comparison  of  rate  constants  determined  from  both  homogeneous  alkaline 
hydrolysis  experiments  and  from  heterogeneous  experiments  (conducted 
with  an  excess  of  S-CL20  in  solution)  indicated  that  rate  constants  were  ap¬ 
proximately  four  times  larger  for  the  homogeneous  reaction  compared  to 
the  heterogeneous  reaction.  Key  here  are  the  findings  that  the  homogeneous 
reaction  rates  were  significantly  larger  than  the  heterogeneous  reaction 
rates  indicating  that  species  with  greater  aqueous  solubility  maybe  better 
suited  to  treatment  by  alkaline  hydrolysis.  The  influence  of  pH  is  also  signif¬ 
icant,  with  higher  pH  leading  to  a  greater  concentration  of  free  hydroxide  in 
the  solution  to  participate  in  the  hydrolysis  reaction.  Thus,  alkaline  hydroly¬ 
sis  may  have  the  potential  for  treatment  of  NTO-contaminated  waste.  It 
may  also  be  effective,  but  to  a  lesser  degree,  for  those  species  exhibiting 
moderate  aqueous  solubility,  including  DNAN,  DNP,  MNA,  MTNI,  and 
TATB.  It  is  essential  for  the  validation  of  the  predictive  approach  to  incorpo¬ 
rate  such  process  variables  as  pH,  or  the  presence  of  inorganic  salt  (salting 
out)  to  adequately  capture  their  impact.  Examination  of  the  susceptibility  to 
nucleophilic  attack  of  particular  bonds  in  the  structure  of  the  energetic  ma¬ 
terial  may  also  provide  valuable  information. 

Chen  et  al.  (2007a)  examined  recovery  of  nitrotoluenes  from  an  aqueous 
stream  using  toluene  as  the  solvent.  They  examined  numerous  process 
variables  including  temperature,  solvent/ wastewater  ratio,  wastewater 
acidity  and  agitation  time.  They  demonstrated  that  recovery  was  strongly 
influenced  by  extraction  temperature,  and  increased  as  temperature  in¬ 
creased.  They  also  found  that  low  pH  and  a  solvent/wastewater  ratio  of  2 
enhanced  recovery.  In  a  subsequent  effort,  Chen  et  al.  (2007b)  also  exam¬ 
ined  the  enhancement  of  the  extraction  process  through  the  addition  of 
inorganic  salts,  and  found  sodium  chloride  to  be  most  effective.  The  use  of 
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solvent  extraction  relies  on  the  affinity  of  a  compound  for  the  organic 
phase;  thus,  extraction  may  be  a  suitable  treatment  method  for  DNAN, 
MNA,  and  S-CL20,  and  to  a  lesser  extent,  DNP  and  MTNI.  However,  with 
the  affinity  of  both  NTO  and  TATB  for  the  aqueous  phase,  solvent  extrac¬ 
tion  would  not  be  suitable  for  remediation  of  streams  contaminated  with 
these  energetic  materials. 

Adsorption  of  energetic  materials  using  activated  carbon  or  similar  mate¬ 
rials  was  also  examined  in  the  literature  (Rodgers  and  Bunce  2001, 
Tomaszewski  et  al.  2003,  Rajagopal  and  Kapoor  2001).  However,  a  review 
of  these  and  other  efforts  indicate  that  assessment  of  adsorption  as  a  re¬ 
mediation  method  for  energetic  materials  is  not  solely  a  function  of  the 
structure  of  the  energetic  material.  Indeed,  the  structure  of  the  adsorbent 
and  the  implementation  of  the  adsorption  process  influence  the  process 
performance. 

Table  10  lists  potential  treatment  technologies  for  waste  streams  contain¬ 
ing  individual  energetic  materials.  The  potential  for  a  particular  treatment 
technology  was  assessed  on  the  basis  of  the  predictions  carried  out  in  this 
work  for  the  various  physicochemical  properties  of  the  energetic  materials. 


Table  10.  Potential  of  existing  treatment  technologies. 


Compound 

Treatment  Approach 

Chemical  Processes 

Physical  Separation  Processes 

Bioremediation 

Oxidation 

Reduction 

Adsorption,  Extraction, 
Membrane 

DNAN 

AH  (M) 

Ex  (H);  Ad(l) 

1 

DNP 

AH  (M) 

Ex  (M);  Ad(l) 

1 

MNA 

AH  (M) 

Ex  (H);  Ad(l) 

1 

MTNI 

AH  (M) 

Ex  (M);  Ad(l) 

1 

NTO 

AH  (H) 

Ex  (L);  Ad(l) 

1 

TATB 

AH  (M) 

Ex  (L);  Ad(l) 

1 

L  -low  possibility;  M-Medium  possibility;  H-High  possibility;  I  -  Insufficient  information 
Ex  -  Extraction;  Ad  -  Adsorption;  AH  -  Alkaline  Hydrolysis 
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5  Conclusion 

This  work  used  the  molecular  structures  of  six  compounds  of  interest  and 
three  reference  compounds  to  estimate  the  physical  properties  including 
critical  temperature  and  pressure,  normal  boiling  point,  vapor  pressure, 
and  enthalpy  of  fusion  that  are  essential  when  using  group  contribution 
methods  for  estimation  of  partition  coefficients.  TURBOMOLE  (a  quan¬ 
tum  chemical  software  package  that  allows  DFT  computations  to  be  readi¬ 
ly  performed)  was  used  to  optimize  geometry  for  each  compound’s  molec¬ 
ular  structure.  The  result  of  the  DFT  computation  was  an  input  file 
containing  a  numerical  description  of  the  electrostatic  potential  surface  of 
the  given  compound. 

This  work  concluded  that,  when  these  electrostatic  potential  surface  de¬ 
scriptions  are  incorporated  into  the  COSMOtherm  program  along  with  es¬ 
timated  physical  property  values,  it  is  possible  to  a  priori  predict  solute 
partitioning  into  the  aqueous  phase  (aqueous  solubility),  between  organic 
and  aqueous  phases  (octanol-water  partition  coefficient),  into  the  gas 
phase  (vapor  pressure)  and  between  gas  and  aqueous  phases  (Henry’s 
constant).  In  the  majority  of  cases  where  experimental  data  were  available, 
the  a  priori  estimates  were  within  experimental  error,  excepting  for  vapor 
pressure  estimates. 

This  work  recommends  that  further  efforts  be  directed  at  experimental  de¬ 
termination  of  vapor  pressure,  aqueous  solubility,  partition  coefficients 
over  a  wider  range  of  temperatures,  so  that  more  meaningful  assessment 
of  waste  treatment  technologies  may  be  realized. 
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Acronyms  and  Abbreviations 


Term 

SDellout 

ACS 

American  Chemical  Society 

AMD 

Advanced  Micro  Devices,  Inc. 

CAS 

Chemical  Abstracts  Service 

CERL 

Construction  Engineering  Research  Laboratory 

COSMO 

Conductor-like  Screening  Model 

DFT 

density  functional  theory 

DNAN 

2.4-dinitroanisole 

DNP 

3,4-Dinitropyrazole 

ERDC 

Engineer  Research  and  Development  Center 

HMX 

octahydro-l,3,5,7-tetranitro-l,3,5,7-tetrazocine 

MNA 

N-methyl-4-nitroaniline 

MSDS 

Material  Safety  Data  Sheet 

MTNI 

l-Methyl-2,4,5-trinitroimidazole 

NTO 

3-Nitro-l,2,4-triazol-5-one 

PDB 

protein  databank 

QSPR/QSAR 

quantitative  structure-property/activity  relationships 

RDX 

cyclotrimethylene-trinitramine 

SMILES 

Simplified  Molecular  Input  Line  Entry  Specification 

SP 

solvation  property 

TATB 

2, 4, 6-Triamino-l,  3, 5-trinitrobenzene 

TNT 

trinitrotoluene 

TR 

Technical  Report 

TZVP 

triple  zeta  valence  Gaussian 

URL 

Universal  Resource  Locator 

WWW 

World  Wide  Web 
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Appendix  A:  Reference  Compound  8-CL20 
*.pdb  and  *.xyz  file  contents. 

CL20.PDB  File  Contents 
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5 
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6 
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7 
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0 
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8 

C 

0 
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9 
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0 
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10 

N 

0 
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11 

O 

0 

-2.545 
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ATOM 

12 

O 

0 
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0+0 

ATOM 

13 

N 

0 
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14 

C 

0 
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15 

N 
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16 

N 

0 
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-1.400 
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17 

O 

0 
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18 
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19 
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O 
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29 
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36 
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1 
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0 
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41 
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2 

1 

3 
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2 

19 

4 

0 
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43 
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4 
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9 
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44 
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4 

20 
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45 
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6 

5 

7 

27 

0 

NONE 

46 
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7 

6 
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47 
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8 

7 

9 

13 
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48 
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9 

8 

4 

10 

0 
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49 
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10 

9 

11 

12 

0 
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50 
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11 

10 

0 

0 

0 
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51 

CONECT 

12 

10 

0 

0 

0 
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52 
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13 

8 

14 

24 

0 
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53 

CONECT 

14 

13 

15 

19 

35 
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54 

CONECT 

15 

14 

7 

16 

0 
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55 

CONECT 

16 

15 

17 

18 

0 
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17 

16 

0 

0 

0 
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57 
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18 

16 

0 

0 

0 
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58 
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36 
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20 

19 

5 

21 

0 

NONE 

60 

CONECT 

21 

20 

22 

23 

0 

NONE 

61 

CONECT 

22 

21 

0 

0 

0 

NONE 

62 

CONECT 

23 
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0 
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67 
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68 
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27 

0 

0 

0 
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69 
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2 
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0 
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70 
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71 

CL20.XYZ  File  Contents 
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Energy : 
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0 

-2.42036 

-2.52717 

1.76569 

N 
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N 
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N 
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0 

-2.54536 

2.88983 

0.96369 

0 
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N 
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-0.70617 
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C 
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N 
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N 
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0 
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0 
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C 
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N 
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-0.53317 
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N 

2 .48164 
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0 

3.34864 

-0.27217 
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0 
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-2.38436 
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N 
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Appendix  B:  Cranium  Predictions  -  Output  for 
Compounds  of  Interest  and  Reference 
Compounds 


HNIW 


Descriptive  Information:  ■  Synonyms: 


Structure  Information:  ■  Molecular  Structure: 


Properly 

Data 

Ref 

Ests 

Ref 

Units 

Mw 

— 

— 

4  3819E-02 

[000] 

— 

dpm: 
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Cm 

rg: 
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m 

Tm: 
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K 

Tb: 

8  6200E+02 

[002] 

8  6176E+02 

[002] 

K 

Tc: 

— 

— 

1  0580E-03 

[003] 

K 

Pc: 

... 

... 

4  8902E+06 

[004] 

Pa 

Vc: 

... 

... 

2.2274 E-03 

[005] 

m3/kg 

Zc: 
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5.4258E-01 

[006] 
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act: 

... 
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2  5472E+00 

[007] 
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— 

— 

— 

— 
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S.298: 

... 
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J/kg 
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— 
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— 

— 

— 
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RDX 


Descriptive  Information: 


Formula:  C3H6N606 
CAS  #:  121-82-4 
Similar  To  -  -  - 


Synonyms: 


Structure  Information: 


Sym  Num: 

Opt  Iso  — 

Chemical  Families: 
Oxygenated 
Unsaturated 


Molecular  Structure: 


N  CH  N. 

I  I 
CV  /CH= 


Constant  Properties: 

Property 

Data 

Ref 

Ests 

Ref 

Units 

Mw: 

— 

— 

2  2212E+02 

[000] 

— 

dpm: 
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... 

... 
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Cm 

rg: 

... 
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m 

Tm: 

— 

— 

— 

— 

K 

Tb: 

6  2660E+02 

[002] 
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K 

Tc: 

— 

— 
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K 

Pc: 

... 
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Pa 
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Vol% 
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Dinitro  pyrazole 
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N-methyl-p-nitroanisole 

Descriptive  Information: 

i 

Synonyms: 

Formula:  C7H8N202 

MNA 

CAS  *  99-52-5 

2-methyl-4-nitro-benzenamine 

Similar  To:  — 

2-amino-5-nitroto!uene 

Structure  Information: 

Molecular  Structure: 

Sym  Num:  -  -  - 

Opt  Iso:  -  -  - 

o.  yo 

V 

I 

r. 

Chemical  Families: 

CH^^CH 

Amine 
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1 

N- 
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1 
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Ref 

Ests 
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Mw 

— 

- 

- 
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— 
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- 
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— 

- 

- 
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K 
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- 
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K 
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— 

. 

. 
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[003] 

K 
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... 

- 

- 
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[004] 

Pa 
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- 

- 
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- 

- 
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[006] 

... 

act: 

... 

- 

- 
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[007] 

... 

Cp,298: 

... 

- 

- 

1  0160E-03 

[008] 

J/kgK 

Hf,298: 

... 

- 

- 

52552E+05 

[009] 

J/kg 

S,298: 

... 

- 

- 

... 

... 

J/kgK 

Gf.298: 

... 

- 

- 

1.6443E-06 

[010] 

J/kg 

Hvb: 

... 

- 

- 

3  7552E+05 

[012] 

J/kg 

Hm: 

— 

- 

- 

1  5707E-05 

[011] 

J/kg 

Hc,298: 

... 

- 

- 

... 

... 

J/kg 

Tf: 

... 

- 

- 

— 

— 

K 

LFL: 

... 

- 

- 

1.2745E-00 

[013] 

Vol% 

UFL: 

... 

• 

• 

— 

... 

Vol% 
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IH-lmidazole,  1  - methyl-2,4, 5-trinitro- 


Formula  C4H3N50© 
CAS#  1©  183-20-1 
Smilar  To:  — 


Imcazole.  1 -methyl-2.4. 5-trnitro- 
1  -Methyl-2.4 . 5-trinitro  midazde 
MTNI 


Structure  Information: 


Molecular  Structure: 


Sym  Num:  — 

Opt  Iso:  — 
Chemical  Families: 
Oxygenated 

Unsaturated 


o 

II 

o=s 


w 

V  / 


o 

II 

N=0 


0=N' 

II 

O 


"N 


CH- 


Constant  Properties: 


Property 

Data 

Ref 

Ests 

Ref 

Units 

Mer 

2.17106+02 

roooi 

... 

cpm: 

... 

— 

Cm 

V 

--- 

— 

m 

Tm: 

7.26026+02 

;ooi] 

K 

Tb: 

6.2S07E+02 

w 

K 

Tc: 

8.450 1E+02 

;oo3] 

K 

Pc: 

5.4707E+06 

;004] 

Pa 

Vc: 

2  2778E-03 

;oo5; 

m3/kg 

2c: 

3.8507E-01 

[006] 

— 

acf 

7.0060E-01 

:oot; 

... 

Cp.238 

... 

— 

JXg  K 

Hf.208: 

... 

— 

J\Q 

S.208: 

... 

... 

J/kgK 

Gf.238 

--- 

--- 

J*0 

Hvb: 

4.24136+05 

;oo8; 

J*0 

Hm: 

— 

— 

Jftg 

Hc.298: 

- -- 

--- 

J\Q 

Tt 

... 

... 

K 

LFL: 

2.W7E+00 

;ooe; 

Vol  % 
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Nitro  Triazolene 


Constant  Properties: 


Property 

Data 

Ref 

Ests 

Ref 

Units 

Mw 

— 

— 

1  3006E-02 

[000] 

— 

dpm: 

--- 

... 

... 

... 

Cm 

rg: 

— 

... 

... 

— 

m 

Tm: 

— 

... 

6  3365E-02 

[001] 

K 

Tb: 

... 

... 

5 684 IE-02 

[002] 

K 

Tc: 

— 

— 

8  2852E-02 

[003] 

K 

Pc: 

... 

... 

8  6052E-06 

[004] 

Pa 

Vc: 

... 

... 

2.2105E-03 

[005] 

m3/kg 

Zc: 

... 

... 

3.5916E-01 

[006] 

... 

acf: 

... 

... 

7.9831  E-01 

[007] 

... 

Cp,298: 

... 

... 

7  6850E-02 

[008] 

J/kgK 

Hf,298: 

... 

... 

1  8991E-04 

[009] 

J/kg 

S,298: 

... 

... 

... 

... 

J/kgK 

Gf.298: 

... 

... 

1  6176E-06 

[010] 

J/kg 

Hvb: 

... 

... 

4  7765E-05 

[012] 

J/kg 

Hm: 

— 

— 

2  7550E-05 

[011] 

J/kg 

Hc.298: 

... 

... 

— 

... 

J/kg 

Tf: 

... 

... 

— 

... 

K 

LFL: 

... 

... 

4  9950E-00 

[013] 

Vol% 

UFL: 

— 

— 

— 

— 

vd  % 

Triaminotrinitrobenzene 


Descriptive  Information: 


Formula:  C6H6N606 
CAS  #:  3058-38-6 
Similar  To:  — 


Synonyms: 


TATB 

1 ,3,5-triamino-2  A6-trinitrobenzer>e 
s-triaminotrimtrobenzene 


Structure  Information: 


Molecular  Structure: 


Sym  Num:  -  -  • 

Opt  Iso:  -  -  - 

°.  yO 

V 

1 

Chemical  Families: 

NHS^c^NH* 

Amine 

II  II 

Aromatic 

//  1  w 

O  NHj  0 

Oxygenated 

Constant  Properties: 

Property 

Data 

Ref 

Ests 

Ref 

Units 

Mw 

— 

— 

2  5815E-02 

[000] 

— 

dpm: 

... 

... 

... 

... 

Cm 

rg: 

... 

... 

... 

... 

m 

Tm: 

... 

... 

9  2651E-02 

[001] 

K 

Tb: 

... 

... 

7.1158E-02 

[002] 

K 

Tc: 

— 

— 

9  1266E-02 

[003] 

K 

Pc: 

... 

... 

56447E-06 

[004] 

Pa 

Vc: 

... 

... 

2.3107E-03 

[005] 

m3/kg 

Zc: 

... 

... 

4.4375E-01 

[006] 

... 

act: 

... 

... 

1  8023E-00 

[007] 

... 

Cp,298: 

— 

— 

94180E-02 

[008] 

J/kgK 

Hf,298: 

... 

... 

31416E-05 

[009] 

J/kg 

S,298: 

... 

... 

... 

... 

J/kgK 

Gf.298: 

... 

... 

1  4887E-06 

[010] 

J/kg 

Hvb: 

... 

... 

4  5020E-05 

[012] 

J/kg 

Hm: 

— 

— 

2  0551E+05 

[011] 

J/kg 

Hc,298: 

... 

... 

... 

... 

J/kg 

Tf: 

... 

... 

... 

— 

K 

LFL: 

... 

... 

1  6717E-00 

[013] 

Vol% 

UFL: 

... 

... 

... 

... 

Vol% 

REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

26-11-2010  Final 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

Prediction  of  Physicochemical  Properties  of  Energetic  Materials  for  Identification  of  Treatment  Technolo¬ 
gies  for  Waste  Streams 

5a.  CONTRACT  NUMBER 

W9132T-06-2-0027 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT 

6.  AUTHOR(S) 

Rebecca  K.  Toghiani,  Hussein  Toghiani,  Veera  M.  Boddu,  and  Stephen  W.  Maloney 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Engineer  Research  and  Development  Center  (ERDC) 

Construction  Engineering  Research  Laboratory  (CERL) 

PO  Box  9005, 

Champaign,  IL  61826-9005 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

ERDC/CERL  TR-10-27 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Armament  Research,  Development  and  Engineering  Center  (ARDEC) 

Energetics  and  Warheads  Division 

Picatinny  Arsenal,  NJ  07806-5000 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

ARDEC 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

This  objective  of  this  work  was  to  predict  the  physicochemical  properties  of  six  energetic  materials.  The  energetic  materials  of  interest  were  2,4- 
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trinitramine  (RDX),  octahydro-l,3,5,7-tetranitro-l,3,5,7-tetrazocine  (HMX),  and  8-CL20. Group  contribution  methods  available  in  the  literature  were 
used  for  predictions.  In  addition,  a  number  of  physical  properties  were  required,  including  critical  temperature  and  pressure,  normal  boiling  point,  melting 
point,  heat  of  fusion,  and  vapor  pressure.  These  were  also  predicted  using  the  chemical  structure  of  the  energetic  compounds.  Where  literature  data  were 
available,  comparisons  with  predictions  were  performed. 

Another  purely  predictive  method,  based  on  chemical  structure,  was  also  employed  for  the  estimation  of  physicochemical  properties.  This  method,  based 
on  the  “Conductor-like  Screening  Model  for  Real  Solvents,”  required  quantum  chemical  calculations  performed  using  TURBOMOLE  to  obtain  the  opti¬ 
mized  electrostatic  potential  surface  of  each  energetic  material.  The  resulting  surface  was  then  used  within  the  COSMOtherm  software  to  predict  solution 
behavior,  including  aqueous  solubility,  octanol-water  partition  coefficient,  Henry's  constant,  and  vapor  pressure. 
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